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Abstract 



The Sun provides an excellent natural laboratory for nonlinear phenomena. 
We use motions of magnetic bright points on the solar surface, at the smallest 
scales yet observed, to study the small scale dynamics of the photospheric 
plasma. The paths of the bright points are analyzed within a continuous 
time random walk framework. Their spatial and temporal scaling suggest 
that the observed motions are the walks of imperfectly correlated tracers on 
a turbulent fluid flow in the lanes between granular convection cells. 



1 



Typeset using REVTgX 



In the outer 30% of its radius the Sun's plasma is convectively unstable, and thermonu- 
clear energy from the core is transported by overturning motion to the visible photosphere. 
The convective motions are most readily seen as the photospheric pattern of bright "gran- 
ules" of hot, rising gas, surrounded by a network of "intergranular lanes" of cooler, darker, 
downflowing gas. The horizontal extent of a granular cell is ~1000 km with a range of 200 
km to 2500 km. The typical separation of granule centroids is 700 km; typical lifetimes are 
5-15 min. 

Magnetic flux threading the solar surface is not uniformly distributed. By processes 
which are not well understood, it is concentrated into discrete kilogauss fibrils of diameter 
< 150 km. These accumulate in the intergranular lanes and vertices where they disturb the 
radiative properties of the photosphere and become apparent to observers through spectral 
line emission. Such magnetic bright points (MBPs) serve as tracers of the evolution of the 
granular flow field and of small scale motions in the intergranular network. 

With Reynolds number « 10 7 turbulent flows are expected in the intergranular lanes. 
Numerical simulations of convection near the solar surface [0 indicate that turbulent vortices 
are formed in the intergranular lanes where the outflows from neighboring granules collide 
and turn downward. Observationally, spectrograms taken in magnetically and temperature 
insensitive spectral lines show line broadening in excess of known Doppler contributions, 
indicating the presence of turbulent flows on the borders of granules 0. These observations 
cannot rule out other types of motion. However, the unique, high-resolution data set de- 
scribed below does provide access to the actual dynamics of the flows and reveals some of 
their quantitative characteristics. 

Here we study subgranular motions at the highest resolution yet achieved The 
data were taken on October 5, 1995 with the 50 cm Swedish Vacuum Solar Telescope on La 
Palma, Spain and consist of two cospatial and simultaneous (within 10 ms) time series of 
178 digital image pairs spanning 71 minutes. One set of images was taken through a filter 
passing wavelengths 430.5 ± 0.5 nm located in the G-band emission of the CH radical. The 
other set of images was made through a wide-band filter passing wavelengths 468.6 ± 5.0 
nm emphasizing the background continuum emission. A 21 Mm x 21 Mm area of enhanced 
magnetic activity near the center of the solar disk was imaged, containing hundreds of MBPs. 
Exposure time was 20 ms, and the pixel scale in the digital images was 60 km on the Sun. 
The images were checked in real time for seeing sharpness, and the best image pair in each 
25 s interval was recorded. Both time series were restored to near the telescope diffraction 
limit (0.2 arcsec) by phase diversity reconstruction ||. The final resolution is 165 km in 
the G-band and 180 km in the broad-band channel. The images were further processed to 
remove instrumental effects and geometric distortion due to atmospheric seeing. 

Next, the wideband images were subtracted from their G-band counterparts, removing 
most of the granulation pattern common to both sets. Finally the images were subjected to 
a high-pass spatial filter to emphasize small-scale features || , and a threshold was imposed 
to identify MBPs. The result is a three dimensional binary (MBP or not) array with two 
350 pixel axes corresponding to x and y spatial directions and a third axis of 178 pixels 
corresponding to the time. A given MBP is tracked by segmentation of its binary object 
tree from the data cube. The segmentation is implemented using one-dimensional voxel 
addressing and searching over limited volumes around a given voxel for valid neighbors. 
After segmentation the individual tree is split into branches and the motion of each sub- 
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object is measured with respect to a common reference point. 

This procedure delivers strings of x and y positions and times, at roughly 24 s intervals, 
for 1800 individual MBPs. The longest sequence for an individual MBP object spans 71 
min, the shortest 9 seconds. During approximately 1 hr the MBPs typically cover an rms 
distance on the order of a granular diameter. Such an MBP path is illustrated in Fig. 
1. The MBPs undergo significant morphological changes on a time scale of 100 s and are 
continuously merging and splitting as a result of their interaction with the granular flow J7| . 
Consequently, the data strings contain a number of large jumps, many representing artifacts 
related to the merging and splitting of two or more MBPs. For example, the trajectory 
shown in Fig. 1 contains four jumps with v > c; the largest is v — 14.6 km s^ 1 . To ensure 
their removal, we follow past practice || and employ a cutoff above the photospheric sound 
speed c = 7 km s _1 . 

Because of the discrete structuring of the magnetic flux, the motions of the MBP tracers 
are naturally described in terms of random walks ||. We will use a general equation for the 
time evolution of the q — th moment of the tracer displacements: 

< r q > 1/q = D ■ r {q)/2 (1) 

where r is the total displacement during time t, and an appropriately weighted average is 
taken over all observed displacements. Fig. 2 shows < r q > 1//,? as a function of t for two 
values of q. We see excellent power law scaling over two orders of magnitude, up to 35 
minutes in time. This permits precise determinations of j(q), which are shown in Fig. 3. 
While Brownian walks give 7(g) = 1 for all q > —1, we find that 7(g) > 1 and varies with 
q. Although diffusion studies are usually limited to mean square displacements, i. e. q = 2, 
here the general q dependence proves to be crucial for correct analysis of the motions. 

Figs. 2 and 3 indicate that for the MBP motions 7(2) = 1.13 ± 0.01 implying superdif- 
fusion, and that in the limit of small q we approach the even more super diffusive value 7(0) 
= 1.27 ±0.01. Our analysis below will indicate that in the ideal situation of infinite, precise 
data we should find that 7(g) is independent of q. The problem here is to find that ideal 
value from limited, imperfect data. When we use no velocity cutoff, a case not illustrated, 
we find 7(2) < 1 indicating subdiffusion, but again that 7(0) ~ 1.3. We will see below that 
it is 7(0), rather than 7(2), that reliably characterizes the MBP walks. 

The super diffusion we have found fits naturally into the framework of the continuous 
time random walk (CTRW) approach ||. This method focuses on the statistical properties 
of the processes under consideration. It handles super diffusion by means of Levy flights, with 
an infinite mean squared displacement, by adding a time cost to the long jumps, so that the 
mean squared displacement as a function of time converges. In a broader context, previous 
experimental work [|IIJ has shown that tracer motion in a rotating flow presents anomalous 
diffusion which can be analyzed in terms of Levy flight motions as we do here. More recent 
experimental work [|TIJ finds that the separations of particle pairs in a turbulent flow, in the 
Richardson regime, show a stretched exponential form rather than the Levy power law form. 
This result, however, can be recovered from particle pairs undergoing relative motions with 
an asymmetric Levy distribution W2 . 



Our use of the CTRW technique is based on two adjustable parameters. The first, A, is 
the exponent of a power law distribution ~ r~ A of spatial steps of tracers due to underlying 
fluid motions. The second, v, connects the spatial step lengths r to elapsed times according 
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to r oc t u . By looking at the spatial scaling properties of the tracer motions, we can directly 
estimate these parameters for our case. Fig. 4 shows the distribution of lengths for all solar 
MBP steps between consecutive images excluding those with v > c. A fit for all r > 40 km 
gives A = 2.65±0.13 ~ 8/3. The assumed step time scaling gives a velocity v = r/t oc r 1-1 ^. 
We estimate this exponent by calculating the velocities of all solar MBP steps with t < 1000 
s, but rejecting those with v > c. A plot of v versus r is shown in Fig. 5. The best fit for 
20 km < r < 1000 km gives 1 — 1/v — 0.35 ± 0.03 ~ 1/3. Thus we find v oc r 1 / 3 , leading 
to v ~ 1.5 ± 0.2 and implying Kolmogorov "K41" scaling with r 2 oc t 3 fl3|Jl4]| . However, we 
have found above that the MBP random walks give < r 2 >oc t 1 with 7 super diffusive but less 
than 3. A key result of the CTRW approach is that, because the tracers may be imperfectly 
correlated with the underlying fluid flow, "Kolmogorov scaling does not necessarily imply 
Richardson's law fl5|." 

To connect our notation with that common in the CTRW literature P JT5| , p!6|| we note 
that the formalism is based on a probability density function ip(r,t) = Cr~^5(r — t u ). Our 
parameter A = — 1/v — (d — 2). In the present case the dimension d = 2. The 1/v term 
arises from the measure associated with the delta function argument. 

The CTRW technique allows analytical calculation of < r q > 1 / £ 'oc t 7 ^)/ 2 for q = 2. The 
time evolution of the moments, that is, the value of 7, is strongly dependent on the values 
of A and v. There are four different cases depending on the degree to which the tracers we 
follow are coupled to the fluid motions. These are illustrated in Fig. 6 as regions in the 
two dimensional A, v plane. In Case I with A > 3 and v > 1/(A — 1) one finds 7(2) = 1. 
This corresponds to a brief correlation time r of the tracers with the flow and thus normal 
Brownian diffusion. In Case II A > 3 but v < 1/(A — 1). This corresponds to strong 
coupling with r = 00 and 7(2) = z/(A — 1) < 1 giving subdiffusion. In Case III A < 3 and 
v > 1/(X — 1). In this intermediate case r is finite but includes significant coupling of the 
tracers with the flow. We find 7(2) = 1 — v(X — 3) > 1 and super diffusion. Finally, in Case 
IV we have 1 < A < 3 and v < 1/(A — 1). Now r = 00, and the tracers are completely 
coupled to the flow. We find 7(2) = 2v. This is super diffusive when v > 1/2. 



The analytical calculation of 7(2) is verified by numerical simulation []16 |. We extend the 



simulations to arbitrary order in the range 0.1 < q < 4.0 by adopting a simple, repetitive 
algorithm for generalized random walks: step in a random direction by a distance r chosen 
with a probability oc r~ A while advancing the time by t oc r 1 ^. 

In the limit of very many trials we find that 7(g) —> 7 = const, which differs from the 
decreasing 7(g) shown by the solar MBPs in Fig. 3. However, the data set is limited in 
time and number of tracers so the statistics are important. Let 7(g) represent the diffu- 
sion exponent calculated by averaging over 100 walks simulated with our CTRW algorithm, 
roughly matching the observational situation. An average over 1000 such 7(g) gives a con- 
stant < 7(g) >~ 1.3, as illustrated in Fig. 3 for a simulation with A = 8/3 chosen to match 
the solar data. Also shown in Fig. 3 are the one standard deviation limits of the 1000 
runs. These indicate the range within which one expects to find a single realization of 7(g). 
Unless gCl, the form of an individual 7(g) is dominated by a small number of the largest 
jumps. Depending on their timing in individual runs they can give large fluctuations in, for 
example, 7(2). For g< 1 the values of 7(g) are determined more equally by all jumps and 
characterize the asymptotic mean of 7(g) for all g. The distribution of the underlying 7(g) 
is narrow for g < 0.5 permitting a robust measurement. 
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The values of A and v determined from the scaling of the solar MBP space and time 
steps are shown on the plot of Fig. 6. These place the solar flows on the v — 3/2 line 
corresponding to Kolmogorov scaling. However, the value of A selects the super diffusive 
intermediate regime between Brownian and Richardson diffusion. Thus, the MBPs are 
imperfect tracers of the flow. For the Richardson case with 7 = 3 we would need perfect 
correlation of the MBPs with the flow. Instead the parameters predict a diffusion exponent 
7 = 1.54±0.39. From the simulations based on the same A and v we find 7(0.1) = 1.28±0.05, 
and 7(2) = 1.24 ± 0.33; the uncertainties show why using q — > is preferable to using 
q = 2 when dealing with data. By direct measurement of the random walks we found 
7(0) = 1.27 ± O.Oland 7(2) = 1.13 ± 0.01. Mindful of the limited data set, we adopt the 
larger value. The result 7 ~ 1.3 holds independent of any velocity cutoff. 

It is no surprise that the MBPs are imperfect tracers of the fluid flows. They are not truly 
local objects, but rather the intersections with the photosphere of extended flux tubes which 
are influenced by various motions at any given time and are restricted to intergranular lanes. 
Nevertheless, the CTRW approach to anomalous diffusion via Levy flights reveals quanti- 
tative properties of the small scale photospheric flows and their coupling to the magnetic 
flux tubes. The spatio-temporal scaling of the MBP dynamics indeed indicates the presence 
of turbulence in intergranular lanes. Our quantifying of the deviation of the MBP motions 
from the underlying fluid motion should contribute to a better physical understanding of 
flux tubes and the solar medium. 
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Berger was supported by NASA SR&T contract NASW-98008. The SVST is operated by the 
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of the Instituto de Astrofisica de Canarias. 
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FIGURES 



FIG. 1. Typical solar MBP trajectory. Distances are in Megameters. Positions were observed 
at intervals of about 24 s; the total elapsed time is 50 min. 

FIG. 2. Moments < r q > 1 / IJ of the cumulative solar MBP displacements r in Megameters versus 
time in minutes for a sonic velocity cutoff and moment orders q = 0.1 (filled circles) and q = 2 
(open circles). If error bars are not visible, they are smaller than the symbols. 

FIG. 3. The multiscaling exponent 7(g) versus moment order q. Solid line: solar MBPs with a 
sonic cutoff v < c, open circles: continuous time random walk (CTRW)mean value for 1000 runs of 
100 walks each, dotted lines: 1 standard deviation limits of the CTRW case indicating the expected 
range for a single 100 walk realization. 

FIG. 4. Filled circles: Logarithms of numbers of step lengths in 5 km bins versus the logarithm 
of step lengths. The solid line has slope -8/3 and represents a fit for all r > 40 km. 

FIG. 5. Filled circles: Step velocities of MBP motions versus step length for all steps occurring 
over 1000 s or less and with v < 7 km s _1 . Error bars are smaller than the symbols. Dashed line: 
best linear fit for step lengths between 20 km and 1000 km. The slope is 0.35 ± 0.03 « 1/3. 

FIG. 6. The two dimensional A, v parameter space. The solid lines A = 3 and v = 1/(A — 3) 
divide the space into four regions with differing correlation properties of the tracers with the 
underlying flows and differing diffusion exponents. The dotted line at v = 3/2 indicates Kolmogorov 
scaling. The heavier dotted line at v = 3/2 and A < 5/3 indicates the locus of Richardson diffusion 
with 7 = 3. The solar symbol represents the values of A and /j, estimated directly from the solar 
MBP walks. 
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